function []=Display(u,Title)
    S = size(u);
    m = min(min(u));
    M = max(max(u));
    step = max(1,floor(S(2)/20));
    clf
%    for i = 1:step:S(2)
%        plot(u(:,i));
%        axis([0,S(1),m,M]);        
%        pause(0.1);
%    end
    plot(u(:,1))
    axis([0,S(1),m,M]);      
    title(Title)
    %set(gca,'nextplot','replacechildren','visible','off')
    set(gca,'nextplot','replacechildren')
    f = getframe;
    [im,map] = rgb2ind(f.cdata,256,'nodither');
    l=length(1:step:S(2));
    im(1,1,1,l) = 0;
    size(im);
    for k = 1:step:S(2)
        plot(u(:,k))        
        axis([0,S(1),m,M]);      
        f = getframe;
        p = (k-1)/step+1;
        im(:,:,1,p) = rgb2ind(f.cdata,map,'nodither');
    end
    name = [Title,'.gif']
    imwrite(im,map,name,'DelayTime',0,'LoopCount',inf)
end